Full Counting Statistics of Avalanche Transport: an Experiment 
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We report the first measurement of high order cumulants of the current fluctuations in an 
avalanche diode run through by a stationary dc current. Such a system is archetypic of devices 
in which transport is governed by a collective mechanism, here charge multiplication by avalanche. 
We have measured the first 5 cumulants of the probability distribution of the current fluctuations. 
We show that the charge multiplication factor is distributed according to a power law that is different 
from that of the usual avalanche below breakdown, when avalanches are well separated. 



PACS numbers: 72.70.+m, 05.40.Ca, 05.60.-k, 85.30.Mn 
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Current noise, i.e. the variance of the current fluctua- 
tions, is the simplest measure of the statistical aspect of 
electronic transport in a conductor, beyond the dc cur- 
rent. Its study as a function of other parameters (voltage, 
temperature, etc.) has been a powerful way to check our 
understanding of the conduction process in many sys- 
tems, and a tool to obtain information that is hidden in 
the mean current 

In order to probe the statistics of the conduction more 
in depth, a better knowledge of the distribution function 
P of the fluctuating current I(t) is necessary, beyond the 
average (/), the dc current, and the variance (i 2 ) with 
i(t) = I(t) — (I). The brackets (.) denote average over 
the distribution P. Most of the time, the full measure- 
ment of P(i) is not possible, but a finite number of the 
moments (i n ) of the distribution can be measured, which 
give some insight into the statistics of i. For example, the 
third moment (i 3 ) reveals the asymmetry (skewness)of 
the distribution around the average. 

Triggered by the calculation of the full counting statis- 
tics of current fluctuations in a quantum conductor Q , 
the measurement of higher moments of current fluctu- 
ations beyond the second has started only recently [!|. 
This measurement and later work 0] have confirmed the 
Poissonian aspect of transport in a tunnel junction, by 
proving that the spectral density of the third moment of 
the current fluctuations is S3 = e 2 (I) (where e is the elec- 
tron charge), independent of temperature. After that, 
measurements of higher moments have been achieved in 
Coulomb blockaded quantum dots @, @. In all these 
measurements, the statistics of transport is driven by the 
finite rate at which electrons can pass a barrier. 

In this article we report the first measurement of 
high order cumulants of current fluctuations in avalanche 
diodes. In such samples, the process that is responsible 
for current fluctuations is not the transport at the one 
electron level, but arises from a collective phenomenon, 
charge multiplication. This occurs due to spontaneous 
creation of electron-hole pairs in semi-conductors in the 
presence of a high electric field, which strongly acceler- 
ates the charge carriers (electrons or holes). When their 



kinetic energy reaches the gap of the semi-conductor, 
they may give part of it to create an new electron-hole 
(e-h) pair which, in turn, will be accelerated and give 
birth to other e-h pairs. Thus, each charge entering the 
sample generates a current pulse of total charge Me with 
M ^S> 1. This avalanche is a statistical process, since 
the e-h pair creation occurs only with a finite probabil- 
ity per unit time: M is not a well defined number but 
has a very wide statistical distribution. This distribution 
has been studied both theoretically and experimentally 
in the regime where no current is injected but individual 
e-h pairs are photo-created, giving rise to well separated 
current pulses @, @, E|. In such a regime, adequate for 
the use of avalanche process as an amplifier of individual 
events, like in the avalanche photodiode, the statistics 
of the pulses determines the noise added by the detec- 
tor. We are not considering here this regime which is 
not stationary, but our goal is to investigate the current 
noise due to avalanche process in the so called breakdown 
regime, when the sample sustains a stationary dc current. 
Semi-conductors working in this regime are indeed used 
as noise sources. 

Experimental setup — The sample we have studied is 
a SM-1 avalanche diode manufactured by Micronetics[loj| 
(we have measured several devices, all showing very sim- 
ilar features). The avalanche starts when the component 
is reverse biased by a voltage > 8V. Too close to this value 
there is no stationary regime. We will focus only on the 
bias region that corresponds to a dc current / > 0.46 
mA where we observe a stationary regime (in the follow- 
ing I = (I) stands for the dc current). The schematics of 
the biasing and measurement circuit is depicted in Fig. 
[lja). We use a bias tee to separate the dc part of the 
current (up to 100 kHz) from its fluctuations (100 kHz 
to 1 GHz) . The sample is biased by applying a dc voltage 
through a i?o = 1 kfi resistor in series with the inductor 
of the bias tee. The fluctuating voltage across the sam- 
ple is detected through the capacitive part of the bias 
tee by a spectrum analyzer (at point A) or amplified (in 
the range 1- 500 MHz), filtered and converted into a dig- 
ital signal (at point D) by a 14 bits, 200MS/s A to D 
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FIG. 1: (a) Schematics of the experimental setup. Zi n is 
the input impedance of the amplifier, that can be modified, 
(b) Spectral density of voltage fluctuations at point A vs. 
frequency for 1 = 1 mA. (c) Equivalent circuit of the sample. 
R s = 4.6 fl, L = 12.6 nH, C = 5.1 pF, R[I = 1mA] = 500 fi. 
i(t) is the current source that models the avalanche. 



converter. All the measurements have been performed at 
room temperature. 

A spectrum analyzer measures the power spectral den- 
sity of the voltage as a function of frequency. This quan- 
tity is proportional to the noise generated by the sample 
at frequency / within a bandwidth of 1Hz. A spectrum 
for a bias of I = 1 mA is shown in Fig. [TJb) . Three 
regions can be distinguished. At low frequency, the noise 
spectral density decays with increasing frequency, ap- 
proximately like 1/f. At higher frequency, the spectrum 
is almost white, up to ~ 500 MHz. In order to remove the 
contribution of the 1/f noise (unrelated to the physics of 
avalanche 11]) and in order to avoid aliasing, we have 
used a 25 — 100 MHz bandpass filter before digitizing the 
signal (indicated by vertical dashed lines in Fig. HUb)). 

In our setup, the voltage vn(t) recorded by the digiti- 
zation card is related to the intrinsic current source i(t) 
that models the avalanche process by vo(f) — H(f)i(f) 
where «£>(/) and i(f) are the Fourier components at fre- 
quency / of vo(t) and i{t). The transfer function H(f) 
contains the effects of the impedances and gains of the ex- 
ternal circuit (bias tee, amplifier, filter) as well as the par- 
asitic impedances inside the component. We have mea- 
sured the frequency response of the setup as well as the 
complex impedance Z s of the sample as a function of fre- 
quency and dc voltage with a vector network analyzer in 
the range 1 — 1000 MHz. We observe that the sample can 
be well modeled by a R s = 4.6 f2 resistor in series with a 
parallel RLC resonator, see Fig. [T{c). R s , L and C are in- 
dependent of the bias voltage, but R varies between 250 il 
and 750 SI when I varies between 0.46 mA and 4.39 mA, 
which implies that H(f) depends on the bias voltage too. 
One also has to know where to incorporate the current 
source i(t) created by the avalanche in the schematics 
of the sample. Knowing that the avalanche itself has a 



FIG. 2: Noise spectral density S 2 (f) = (|i(/)| 2 ), vs. dc cur- 
rent (symbols: measurement, solid line: fit 7 -1 ' 21 ). Inset: 
measured probability distribution of vd for 1 = 1 mA. 



white spectrum Q, we conclude from the frequency de- 
pendence of the noise spectral density (Fig. [IJb)) that 
i(t) must be in parallel with R s , see Fig. [He). The high 
frequency roll-off observed on the power spectrum comes 
from the RLC circuit (l/(2^^^/LC) = 600 MHz). Thus, a 
direct determination of i(t) could be obtained by divid- 
ing by R s the voltage across the diode, measured with 
a high impedance voltmeter. We have used an amplifier 
with 50 f2 input impedance and calculated H(f) (which 
includes the finite input impedance of the amplifier) to 
deduce i(f) from vo(f)- 

Second moment - We have recorded the voltage fluc- 
tuations vo(t) and made histograms of them for many 
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bias points, see inset of Fig. O The variance 
(here the brackets mean time averaging) is related to 
the noise spectral density 82(f) — (\i(f)\ 2 ) by: (uf>) = 
J-™\H(f)\ 2 S 2 (f)df. Since we observe that S 2 (f) is 
reasonably frequency independent within the detection 
bandwidth (see Fig. QJb)), the variance of vd and 
the spectral density of the current fluctuations 5*2 are 
simply proportional: (i>fj) = H 2 S 2 where H 2 (V) — 
J\H(VJ)\ 2 df. We show on Fig. rj the variations of 
S 2 with the dc current running through the sample. 
Note that the amount of noise emitted by the avalanche 
diode is huge, so that the contribution of the amplifier 
is negligible. For I = 1mA, S 2 = 1.15 x 10~ 15 A 2 /Hz 
which is equivalent to the Johnson noise of a 50 fl resis- 
tor at a temperature of 18 millions Kelvin degrees, to 
compare with the 300 K noise temperature of the ampli- 
fier. Defining the Fano factor F2 by S2 = i^eJ, we find 
F2 = 6.8 x 10 6 for I — 1 mA. The noise is thus extremely 
super-Poissonian (F2 = 1 for Poisson noise). 

Third moment and environmental effects — From the 
histograms of vo(t) we have calculated the third mo- 
ment (wfj). It is related to spectral density ^(/i,^) = 
(i( fi)i(f2)i(— fi — /b)) of the third moment of i(t) by 
--JjH(f 1 )H(f 2 )H(-f 1 - f2)S 3 (f 1 ,f 2 )df 1 df 2 . 
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FIG. 3: Fano factors F„ = CVi/(e n-1 J) of the first cumulants 
of the distribution of current fluctuations, as the function of 
the second Fano factor F2 = Cn/iel), in log-log scale. 

If we suppose that like S2, S3 is frequency inde- 
pendent [13j, one has: (v^) — H3S3 with H 3 — 
J $dhdhH{h)H(f 2 )H{-h - f 2 ). As it has been 
demonstrated experimentally and well understood [3j, [l2j, 
0| , there is an additional contribution to S3 given by: 
ReffS2^jy- Reff is an effective resistor that contains the 
impedance of the sample and that of the environment, 
here the internal RLC circuit and the input impedance 
of the amplifier Zi n . In order to estimate the environ- 
mental contribution, we have modified R e ff by adding 
resistances to ground at the input of the amplifier, thus 
varying Zi n in the range 4 — 50 SI. This had no effect 
on S3 (besides an overall reduction due to voltage divi- 
sion by the added resistances) except at very low cur- 
rent, where both 5*2 and \dS2/dV\ are maximum, and 
where we observe a 5% negative deviation, in agreement 
with the prediction. Thus, environmental contributions 
to S3 are negligible, i.e. we have access to the intrin- 
sic third moment of the current fluctuations. Moreover, 
5*3 (and higher moments) not being influenced by Z{ n 
demonstrates that the statistics of the current is not af- 
fected by the non-linearity of the component (which is 
small since Lu) <C R in our detection range). We have 
plotted the result for S 3 (I) in Fig. [3l as F 3 = S 3 /{e 2 I) 
vs. F 2 — S2/(eI). We discuss the interpretation of this 
result later. 

Higher order cumulants — Any statistical measure- 
ment has to face a fundamental problem: the central limit 
theorem, according to which the distribution of the sum 
of N independent variables tends to a Gaussian when N 
goes to infinity, regardless of the distribution law of each 
of the variables. In our measurement, a current 1=1 mA 
integrated during 5 ns corresponds to N « 3 x 10 7 elec- 
trons, and the measured histograms are very close to 
gaussians, see inset of Fig. Nevertheless one can access 
the probability distribution P(i) by working with its cu- 
mulants, defined as follows. We introduce the character- 
istic function x( z ) = ( e * z ) (* is the fluctuating current). 



The Taylor expansion of x( z ) — Enl'") 2 ™/' 1 ' provides 
the moments (i n ). The cumulants ((i n )) are defined by 
the Taylor expansion of In x{ z ) = K P(i) 

is Gaussian, lnx(z) is a second degree polynomial, and 
((«")) = for n > 3. Thus the cumulants reveal how 
much a distribution deviates from a Gaussian. From the 
histograms of vrj(t) we calculate the cumulants ((i>jr>)) 
from which we deduce the spectral cumulants of P(i), 
which we note C n , supposing that they are frequency 
independent, as C2 = 5*2 and C3 = S3 are. A simple 
calculation gives indeed: ((«£>)) = H n C n where the coef- 
ficients H n are given by: H n = f^ 00 H n (t)dt with H(t) 
the inverse Fourier transform of H(f), i.e. the impulse 
response of the circuit. This definition coincides with 
that given before for n = 2 and n — 3. In order to check 
the calibration coefficients H n , we have measured the 5 
first cumulants of a pseudo-random current generated by 
a lGS/s arbitrary waveform generator with a statistics 
that we have computed. The ratio of the measured cu- 
mulants to the real ones provides us with another set 
of calibration coefficients H' n . We obtain H n ~ H' n up 
to n = 5, see inset of Fig. HJ The systematic error on 
the sixth spectral cumulant is only by a factor 2. We 
conclude that we can perform a reliable and quantitative 
measurement of C n up to n = 5. 

We have measured the 5 first spectral cumulants of the 
current noise of the avalanche diode. We observe that 
they all decrease as a function of the average current as 
a non- integer power law (see e.g. S2 on Fig. In order 
to synthesize our results, we show on Fig. [3] the Fano 
factors F n of the n th cumulant (n > 3) as a function of 
the second one F2; we define the (dimensionless) Fano 
factors by: C n = F„e" -1 /. We observe a remarkable 
behavior: F n {I) oc F^ 2 (I), which can be rewritten as 
F n (I) = g„G n (I) (n > 2) where the function G(I) is the 
same for all the cumulants and g n are constants. Note 
that all the cumulants can be in principle calculated from 
the data, but the uncertainty on C n due to the statistical 
aspect of the measurement (we record i>i)(i) during a 
finite time) grows fast with n [lj| . Fig. [3] corresponds 
to 8 x 10 10 independent measurements of the current for 
each bias point. 

Interpretation - The result of Fig. [3] is central in 
the understanding of the statistical aspect of transport 
in the avalanche regime. A Poisson distribution with an 
effective charge q = Me, has Fano factors F n = M n ~ x = 
F^ 1 , even if the effective charge depends on /. Our mea- 
surements clearly rules out this scenario. The reason for 
this is that the avalanche itself is responsible for the cur- 
rent noise. We will now compare our experimental result 
with what is predicted in the non-stationnary regime. An 
injected charge gives rise to m charges with a probabil- 
ity P(m). Thus, an injected photo-current io gives rise 
to a dc current I = (m)i and to current fluctuations of 
spectral cumulants C n = e n_1 io((m p )) [1, H^. The mul- 
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FIG. 4: Ratio of Fano factors F3/F2 = Cz/(eC2) as a func- 
tion of the dc current I, in log-log scale. Inset: ratio of two 
independent sets of calibration coefficients for the cumulants. 

tiplication factor to has been shown to obey a probability 
distribution such that: ((to")) ~ (to™) cx (to) 2 ™" 1 [IGif. 
This implies for the Fano factors: F n cx (to) 2 ™ -2 cx 
in strong contrast with our results. 

The non-stationnary regime is a particular case of 
avalanches, which are more generally described by a 
probability distribution that is a power law, P(m) ~ m~ T 
(with 1 < t < 2) valid for to 1 ~ 1 < to < m-> fl7l]. 
The results of the non-stationnary regime are recovered 
taking r = 3/2 [8], the mean field result of theory of 
avalanches . toi represents the most probable avalanche, 
which is small, and m.2 the largest one, which is some- 
what rare but dominates the statistics of m. Indeed, 
((m p )} ~ (m p ) cx m.2 +1-T . The average multiplication 
factor is (to) ~ to 2 , -1 " -C mi, and the current spectral 
cumulants C p — k p e p ~ 1 iom P , +1 ~ T . The coefficients k p 
depend on the exact shape of P(m) for to < mi and 
to > TO2 . This describes our experimental results pro- 
vided I cx ioTO,2" r . Then the function G(I) is sim- 
ply G(I) = TO2- However, in the stationary regime, in 
which our experiments have been performed, the dc cur- 
rent is imposed by a battery, and not related to any io, 
which should disappear from the theory. This implies 
a power law relation between the average current and 
the largest avalanche m^- I cx m 2 T , which reflects how 
the probability distribution of the multiplication factor 
adapts to the applied dc current. From our experimen- 
tal data we can extract the current dependence of mi 
through F3/F2 — {kz/ki)mi. As can be seen in Fig. |4j 
we indeed observe a power law, from which we deduce 
t ~ 1.87. The largest avalanche is huge, mi = 1.38 x 10 6 
for 1=1 mA and ki — k% , whereas the average avalanche 
is modest, (to) = 6.3. From F n = k n m^ and taking 
k-2 = &3 we deduce ki = 3.6 x 10~ 6 , k^ = 8.6ki and 
&5 = 20.9/C2- A detailed theory is of course needed, in 
particular to explain how the avalanche feeds itself to 
reach the stationary regime. 



In conclusion, we have demonstrated how the statisti- 
cal physics of a complex phenomenon can be extracted 
from the analysis of its fluctuations beyond the second 
moment. In particular, we have demonstrated that the 
charge multiplication factor of avalanche diodes has a 
power law probability distribution different from what 
has been reported until now in the non-stationary regime, 
and different from what is predicted using mean field ap- 
proximation in avalanche theories. This may open new 
routes for the study of other components, or complex con- 
ductors such as those with charge density waves. More- 
over, avalanche diodes are advertised as perfect gaussian 
noise sources. This is definitely incorrect and may have 
practical consequences. For example, using the sign of 
the instantaneous fluctuating current as a source of ran- 
dom bits gives, for the sample we have analyzed and in 
our bandwidth, P(i > 0) ~ 1.01P(i < 0). To cure this, 
the central limit theorem may help, with the price of a 
smaller bit rate. 
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